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Objects  and  goals  of  developing  an  advanced  SOG 


This  one-year  EOARD  grant  was  awarded  to  promote  the  co-operation  of  COIL  Laboratories 
at  Institutes  of  Physics  in  Prague,  Czech  Republic  and  Samara,  Russia  on  development  new  concepts 
of  advanced  chemical  generator  of  singlet  oxygen  for  a  COIL.  The  generator  should  have  the 
following  characteristics: 

1.  A  high-pressure  operation  with  ~100  Torr  of  O2  partial  pressure 

2.  A  single  BHP  pass  burn  down  (probably  requiring  a  dilute  BHP  to  maintain  the  thermal  control 
and  so  to  limit  KBO  vapor) 

3.  A  close  coupling  of  the  generator  to  the  nozzle  plenum  to  minimize  CB^A)  losses 

4.  An  efficient  disengagement  of  gas/liquid  (conceptually  involving  the  angular  acceleration  of 
fluid  and  gas  or  both) 

5.  A  scalability  to  be  suitable  for  the  airborne  and  mobile  applications 

The  investigation  within  the  framework  of  this  grant  has  been  aimed  at  conceptualisation  of  the 
problem  including  a  literature  study,  theoretical  estimation  of  particular  problems,  calculations,  and 
designing  the  concept  of  the  experimental  device. 
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A.  The  conceptual  designing  of  the  small-scale 
Rotating  Sparger  SOG  (RSSOG)  (30  mmole/s) 
with  next  requirements 


1.  Estimations  of  hydrodynamic  and  mass-transfer  rates 
chlorine  utilization  and  C^A)  yield 

1.1.  Basic  data  and  formulas 

The  estimation  of  RSSOG  parameters  is  based  on  the  simplified  liquid-gas  hydrodynamic  and 
simplified  model  of  the  transport  phenomena  of  the  bubbler  centrifugal  apparatus  performed  in  the 
interim  Report  00 1  of  this  grant. 

The  parameters  of  RSSOG  that  should  be  found  for  achieving  the  required  parameters  of  the 
advanced  SOG  are  presented  in  Table  1. 

Tab.  1 


1 

R 

Radius  of  RSSOG 

2 

do 

Diameter  of  holes 

3 

H 

Height  of  liquid  layer 

4 

w 

Axial  BHP  velocity 

5 

n 

Revolutions  per  minute  (RPM) 

6 

Ug 

Gas  velocity  from  the  nozzle 

7 

Nh 

BHP  molarity  (KOH) 

Tab.  2 

Physical-chemical  data  of  the  gas  and  BHP 


1 

Pi 

BHP  density 

1 .3  g/cnf* 

2 

Pi 

BHP  viscosity 

0.1  g/cm/s 

3 

K 

Rate  constant  of  reaction  Cf+HOj  ->  02(‘A)+2Cr+H+ 

~5x108M-1s-1 

4 

Hc 

Henry  constant 

~1 

5 

Dc 

Chlorine  diffusion  coefficient  in  BHP 

5xl0"6cm2/s 

6 

Do 

Oxygen  diffusion  coefficient  in  BHP 

5xl0"6cm2/s 

7 

Dh 

HO  2  diffusion  coefficient  in  BHP 

5xl0"6cm2/s 

8 

m 

Chlorine  viscosity 

1.6xl0"4  g/cm/s) 

9 

DcHe 

Chlorine  diffusion  coefficient  in  He 

300/P(torr)  (cm2/s) 

10 

DC02 

Chlorine  diffusion  coefficient  in  O2 

100/P(torr)  (cm2/s) 

11 

Ta 

Lifetime  of  02(!A)  in  the  BHP 

2x10"6s 

5 


Tab.  3 

Given  and  calculated  intermediate  parameters 


1 

T=253°K 

Input  BHP  temperature 

2 

Po=100  torr 

Total  pressure  (02  +  Cl2resid)  over  BHP  layer 

3 

A=0 

Dilution  with  Helium 

4 

wr=t»R 

Peripheral  BHP  velocity 

5 

G=(b2R 

Centrifugal  acceleration 

6 

P=  P0  +co2RHpl 

Gas  pressure  at  the  nozzle  outlet  (torr) 

7 

pB=(P/RT)(71/(l+A)+A*4/(  1+A)) 

Gas  density  at  the  nozzle  outlet  (g/cm3) 

8 

NC=P/RT/(1+A) 

Chlorine  concentration  at  the  nozzle  outlet  (mole/cm3) 

9 

ps=  3x10'4xA/(1+A)+  1.6x10"4/(1+A)) 

Gas  viscosity  (g/cm/s) 

10 

Dc8=( ( 1 +A )/P ) ( A/  300+  1/100)1 

Chlorine  diffusion  coefficient  in  02-He  mixture  (cm2/s) 

11 

Re-  P'U°d|' 

Reynolds  number  for  the  gas  in  nozzle 

Tab.  4 

Estimation  of  the  hydrodynamic  parameters  of  the  rotary  bubbler,  specific  chlorine  absorption  time 

and  02(‘A)  yield  parameters 


1 

t  r..06.  (p.+.c  fun‘“ 

■"  068i+K/„,r[o.oJ 

Radius  of  bubble 

2 

U  » 1.43  rG 

Velocity  of  bubble 

3 

T  '  0. 

c 

C=  ® 

A/dol 

12  x  10~ 

Pgd0G 

.  P^Us  . 

'6C2  + 1.05  x  10'3C  +  0.21  -  0.35  x  10’3 

\°‘ 1 

Re006 

/ 

Time  of  bubble  formation 

4 

f  = 

0.25  x  C2  +  0.71  x  C  +  0.19  x  106  -  0.44  x  103 

Frequency  of  bubbles  generation 

5 

tr  =  x  +  H/U 

Bubble  residence  time  in  BHP 

6 

tg  «0.06r/D8 

Chlorine  absorption  time  due  to  gas 
resistance 

7 

t.«Hcr/( 

'  knhdc 

) 

Chlorine  absorption  time  according  to 
surface  reaction 

8 

tc  =  2r/U 

The  time  of  the  BHP  surface  renewal 

9 

^dep  ~ 

(n„  Y  H2  jxDh 

Time  of  strong  H02  depletion 
(formation  of  sheet  flame) 

DcKNh  4 

10 

t  =  +  [ 
*  N,1 

'  r3  Y 

udhJ 

“( 

o _ 

Time  of  total  chlorine  absorption  in 
sheet  flame  approximation 

11 

Ut  =  -h 
tabu 

1 - 1 

^  |  <n] 

i 

1A) 

Estimated  chlorine  utilization  for 
moving  bubble  (see  below) 

6 


12 

v  nQ  f  2N  7iD  r  ) 

Yjet  ~  0.9exp  -  c  c 

1  HcNh  UD0ta  ) 

Estimation  of  minimal  detachment 

02(1A)  yield 

13 

Pctr=(P+Po)/(A+l)*tr/4 

Average  Px  parameter  due  to 
homogeneous  quenching  of  CE^A) 

14 

Y=Ydet/(l+Yde,Pctr) 

Estimation  of  CBC'A)  yield 

15 

N„=(2r)-2 

Maximum  density  of  chlorine  nozzles 
(1/cm2)  to  avoid  bubble  coalescence 

16 

.  .  ™ .  47ir3  xr 

Mc=fNcx  N„ 

Estimated  oxygen  flux  from  the  bubble 
layer  (mmole/cm2/s) 

Note  to  Item  11:  Estimation  of  chlorine  utilization: 

Because  of  very  high  initial  chlorine  pressure  (more  than  100  torr  under  BHP  layer)  and  very 


fast  depletion  of  the  HO  2  ions,  the  chlorine  absorption  rate  will  occur  according  to  the  flame  sheet 
approximation.  In  this  approximation  the  chlorine  flow  rate  to  bubble  surface  will  be 


J,  =  71— N„D, 


3Ur 

27iD, 


The  chlorine  absorption  from  the  moving  bubble  is  described  by  the  equation 


We  will  assume  that  pressure  in  the  bubble  approximately  equals  to  pressure  in  an  ambient  BHP  layer. 
Hence  the  radius  of  the  bubble  increases  according  to 

r(x)3  =  P/P(x)r3 .  (2) 

where  P=Po+piGh,  P(x)=  P=Po+PiG(h-x).  The  residual  chlorine  partial  concentration  in  the  bubble  is 
given  by 

Nc(x)  =  ^)qc(x)  (3) 

where  qc  is  the  molar  fraction  of  the  residual  chlorine.  The  substitution  of  (2)  and  (3)  gives  equation 
for  r|c 
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NH| 

'  3Dh  ' 

0, 

Po 

Nct 

v2jir3Lb 

J 

_P0  +  pG(h  -  x) 

dilc 

dx 


The  solution  of  this  equation  gives  the  estimation  of  the  chlorine  utilization 


(4) 


Ut  =  2N^3DH 


N„ l 2nr  U 


-1 


P-P„ 


tabu 


-1 


P-Pn 


(5) 


In  this  solution  the  bubble  velocity  was  assumed  constant. 

Indeed  U=1.43(r(x)G)°'5=1.43(rG)05(P/P(x))1/6  and  weakly  depends  on  x.  The  formula  (5)  is  valid 
when  calculated  Ut<l.  The  case  Ut»l  means  that  the  flame  sheet  approximation  is  not  valid  and  a 
surface  reaction  approximation  have  to  be  used.  When  Ut>  1  h-5  .  the  absorption  mode  is  intermediate 
between  flame  sheet  approximation  and  surface  reaction  approximation.  In  this  case  it  is  necessary  to 
compare  values  of  times  ts,  tr. 

Note  to  Item  12:  The  estimation  of  the  detachment  yield  in  the  flame  sheet  model  approximation  is 


Y, 


exp 


f  2Nc  inDp 
V  HcNh  l>c-w 


a  + 


D„ 


i  0.9 


where  t  is  a  time  of  starting  of  the  BHP-Ch  contact.  The  most  of  the  oxygen  gas  is  produced  near  the 
bubble  equator  (t=r/U).  Hence  the  estimation  of  the  Oof' A)  detachment  yield  is 


>  0.9exp 


2N^  7tDcr 

V  hcnh  udo^ 


Note  to  Item  13:  Pi  is  defined  as  an  average  for  pressure  (P  +  Po)  and  divided  by  2  because  the 
CbC'A)  is  monotonically  produced  during  time  tr. 

Note  to  Item  14:  It  is  assumed  that  distance  between  nozzles  should  be  more  than  bubble  diameter. 


1.2.  Estimation  of  RSSOG  parameters  for  achieving  Po=100  torr. 

A  one  example  of  the  estimations  of  RSSOG  parameters  have  been  performed  for  the  next  a 
priory  given  input  parameters  and  varied  number  of  the  revolutions  per  second: 


Tab.  5 


1 

R=3  cm 

Radius  of  RSSOG 

2 

do=0.02  cm=200  pm 

Diameter  of  holes 

3 

H=0.3  cm 

Height  of  BHP  layer 

4 

w=0 

Axial  BHP  velocity 

5 

n-varied 

Revolutions  per  minute  (RPM) 

6 

Uo=130  m/s 

Gas  velocity  from  nozzle 

7 

Nh=6  mole/L 

BHP  molarity  (KOH) 

The  results  of  calculations  according  to  the  formula  in  Table  4  and  for  the  input  parameters  in  Table  5 
are  presented  in  Figs.  1-5.  The  dependences  presented  in  Figures  1-5  can  be  explained  qualitatively. 
When  the  centrifugal  acceleration  G  is  not  too  big  the  pressure  P  is  close  to  Po,  and  according  to  eq. 
(5),  the  utilization  increases  with  increasing  RPS.  At  higher  G,  the  pressure  P  is  substantially  higher 
than  Po  and  the  chlorine  concentration  in  the  nozzle  increases  with  increasing  G.  It  explains  non¬ 
monotonic  dependence  of  the  utilization,  detachment  yield  and  Pi  parameter  with  increasing  RPS. 


Fig.  1  The  initial  bubble  radius  as  a  function  of  the  rotation  rate  (RPS) 
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Fig.  4  Dependence  of  detachment  OiI'A)  yield  (lower  limit)  on  RPS 

A  lot  of  calculations  have  been  performed.  How  the  parameters  indicated  in  Figs.  1-5  are 
changed  when  other  parameters  in  Table  5  are  changed: 

a)  Nozzle  diameter.  A  decrease  in  the  nozzle  diameter  (lower  than  200pm)  results  in  a  smaller  bubble 
size,  lower  bubble  velocity  and  higher  utilization.  However,  simultaneously  Pt  increases  and  Ydet 
drops.  It  is  necessary  to  decrease  simultaneously  the  BHP  layer  height  to  maintain  a  high  utilization 
and  02(*A)  yield. 

b)  The  gas  velocity  in  the  nozzles.  An  increase  in  the  gas  velocity  results  in  increasing  the  bubble 
size  and  bubble  velocity,  which  causes  the  utilization  decrease.  The  maximum  density  of  the  nozzles 
in  the  bubbler  decreases,  the  frequency  of  the  bubble  production  decreases,  and  consequently  the 
maximum  oxygen  flux  doesn’t  increase  substantially  with  increasing  Ug. 

c)  The  height  of  the  BHP  layer.  With  an  increase  in  H,  the  utilization  increases.  However,  increasing 
H  results  in  the  increase  in  the  flight  time  tf,  and  the  increase  in  the  bottom  pressure  P.  This  causes  the 
Ydet  drop  and  Pi  increase. 

d)  The  BHP  molarity.  The  decrease  in  NH  substantially  reduces  the  chlorine  utilization  and  the 
detachment  yield  at  a  high  initial  chlorine  pressure.  To  attain  the  required  parameters  with  diluted 
BHP,  it  is  necessary  to  decrease  the  nozzle  diameter,  gas  velocity  in  the  nozzle,  the  BHP  height,  and  to 
increase  the  density  of  the  nozzles. 
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Fig.  5  Dependence  of  the  quenching  parameter  Pi  on  the  RPS 


Several  sets  of  the  parameters  presented  in  Table  1  have  been  used  in  calculations.  We  select 
only  those  sets  of  the  initial  parameters,  at  which  the  required  output  parameters  have  been  achieved. 
Table  6  presents  the  estimation  of  other  sets  of  the  input  parameters  resulting  in  achieving  a  high 
utilization  and  Oil' A)  yield  of  ~70%  at  the  output  pressure  Po=100  torr. 

Tab.  6 

The  recommended  initial  RSSOG  parameters  for  achievement  of  the  high  utilization  and 
02(*A)  yield;  the  radius  of  the  bubbler  apparatus  is  R=3  cm 


Ini 

put  parameters 

Output  parameters 

n 

s4 

H 

mm 

do 

jam 

Ug 

m/s 

Nh 

mole/L 

Wr 

m/s 

r 

mm 

U 

cm/s 

Ut 

Yde, 

Pi 

torrxs 

Y 

Mc 

mmole/cm2/s 

1 

60 

3 

200 

130 

6 

11.1 

0.8 

250 

~1 

-80% 

0.14 

-0.7 

-1 

2 

100 

1.5 

100 

130 

6 

18.8 

0.4 

300 

~1 

-85% 

0.08 

-0.7 

-1 

3 

80 

6 

300 

130 

6 

9.4 

0.9 

3.8 

~1 

-70% 

0.32 

-0.5 

-3 

4 

100 

1 

50 

50 

2 

18.8 

0.15 

200 

~1 

-78% 

0.06 

-0.75 

-0.5 

5 

80 

1 

50 

100 

3 

9.4 

0.2 

200 

~1 

-80% 

0.06 

0.76 

-0.5 

Note:  For  all  initial  parameters  tg<tab  and  tdep«ts.  Hence  the  flame  sheet  model  approximation  is  valid. 

The  parameters  in  Table  6  correspond  to  only  very  rough  estimation  of  RSSOG  parameters.  The  main 
parameter  for  the  COIL  is  the  chlorine  flux  Mc  or  oxygen  molar  flow  rate  produced  from  1  cm2  of  the 
bubbler  layer.  Parameterization  showed  that  a  high  Ydet>80%  and  low  Pi<0.1  torrxs  could  be  achieved 
at  MC<1  mmole/cm2/s,  which  is  not  enough  large  for  the  use  in  a  COIL.  The  main  reason  is  in 
requirement  of  the  high  pressure  Po=100  torr.  The  quenching  of  CbC'A)  and  depletion  of  H02  in  the 
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bubbler  layer  is  very  strong  to  ensure  simultaneously  a  high  utilization  and  Y>50%  at  MC>1 
mmole/cm2/s  and  Po=100  torr.  Moreover,  it  is  quite  difficult  to  ensure  good  output  parameters  at  the 
low  BHP  molarity.  A  higher  Mc>3  mmole/cm2/s  is  possible  to  achieve  at  lower  predicted  yield 
Ydet<70%  and  higher  Pt>0.2  torrxs.  So  at  higher  Mc  the  predicted  ChC’A)  yield  would  be  lower  than 
50%.  It  seems  that  the  real  liquid-gas  hydrodynamic  and  mass  transfer  phenomena  are  much  more 
complex.  The  heating  of  the  BHP  results  in  increasing  in  the  mass  transfer  coefficients  (Dc,  Dh,  D0), 
the  decrease  in  the  BHP  viscosity.  These  phenomena  have  to  improve  chlorine  utilization.  The  real 
output  parameters  of  RSSOG  can  be  much  better  at  higher  Mc.  But  only  experiment  can  clear 
these  problems  and  questions. 

1.3.  Parameters  of  RSSOG  at  Po<100  torr. 

It  was  found  that  a  decrease  in  the  output  pressure  P0  makes  possible  to  increase  the  RSSOG 
specific  output  flux  Mc  and  to  save  a  high  chlorine  utilization  and  OgC'A)  yield.  Two  sets  of  the  input 
RSSOG  parameters  at  P0=20  torr  are  presented  in  Table  6a. 


Tab.  6a 

Example  of  the  input-output  parameters  of  the  RSSOG  at  moderate  Po=20  torr 


Input  parameters 

Output  parameters 

n 

H 

do 

ug 

Nh 

r 

U 

Ut 

Yde, 

Pi 

Mc 

s’1 

mm 

jam 

m/s 

mole/L 

mm 

cm/s 

torrxs 

mmole/cm7s 

6 

80 

5 

500 

200 

6 

1.6 

500 

~1 

>75% 

0.16 

~3 

7 

80 

4 

300 

100 

3 

0.1 

374 

~1 

70% 

0.12 

~2 

2.  Designing  of  parameters  for  the  single  pass  burn  down  small-scale 
RSSOG  (M(02)=  30  mmole/s,  P0=100  torr) 

Additional  specifications: 

r|w  -  water  vapour  fraction 

PWSi  -  saturated  water  vapour  pressure  for  the  input  BHP  temperature  T 
PWSf  -  saturated  water  vapour  pressure  for  the  output  BHP  temperature  Tf 
Y  -  expected  OdC'A)  yield 
Z  -  axial  length  of  the  rotary  bubbler 

Nuf  -  output  concentration  of  the  HO  2  ions 
tBHP  -  BHP  residence  time  in  RSSOG 
M(02)-  output  oxygen  molar  flow  rate 
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In  the  single  pass  burn  down  RSSOG  the  output  water  vapour  fraction  should  be  enough  low  and  the 
BHP  consumption  high.  The  estimations  made  below  are  based  on  the  results  present  in  the  Report 
001.  We  assume  that  the  water  fraction  r|w=15%  to  the  oxygen  concentration  does  not  decrease 
substantially  the  COIL  gain  and  power  in  comparison  with  r|w=0.  The  fraction  of  the  water  vapour  is 
given  by  ((24),  Report  001): 


(1  +  A)Pwsi  (Pwsf/Pwsi)-1 


In 


V  Pwsi  J 


(6) 


The  water  vapour  partial  pressure  at  input  BHP  temperature  is  given  by 

PwSi0orr)  =  exp^23.44 -  ^ '  j 


(7) 


The  water  vapour  partial  pressure  at  output  BHP  temperature  is  given  by 

Pwsf(t0rr)  =  eXP^23-44-^1j  (8) 

Assuming  T=253K,  dilution  A=0,  and  Po —  1 00  torr,  we  obtain  the  required  output  BHP  temperature 
Tf<315K  to  maintain  qw<0.15.  The  required  maximum  BHP  residence  time  (tenp)  in  the  RSSOG  to 
ensure  Tf<3 15K  and  qw<().  1 5  is  given  by 


t  =  —  -  (t  -  T  ) _ ^bhp _ 

BHP  W  lf  1  Mc(Qr  +  (1  -  Y)Qa) 


(9) 


At  maximum  residence  time  (9)  the  output  H02  concentration  is  given  by 

NHf=NH-2|^tBHP  (10) 

In  (9)  and  (10),  Cbhp=3.7  J/cm3/K  is  the  volume  heat  capacity  of  the  BHP,  Qr=112xl03  J/mole  is  the 
heat  release  in  reaction  CI2+2HO2"— ^CT+HiO+Oat'A),  QA=94xl03  J/mole  is  the  heat  release  in 
reaction  02('A)— >02(3S). 

For  the  oxygen  molar  flow  rate  M(02)  and  radius  of  the  bubbler  cylinder  R=3  cm  the,  length  Z  of  the 
perforated  bubbler  cylinder  in  the  axial  direction,  the  BHP  axial  velocity  w  and  minimal  volumetric 
rate  Vbhp  are  correspondently  equal  to 


z_  M(Q2) 

Mc2jiR 

Z  M(0,) 

w  = - = - - - 

^BHP  2ttRt  bhp 


VBHP  =  2jiRHw  =  H 


M(02)  _  M(02)(Qr  +  (1  -  Y)Qa) 


^c^BHP 


CBHp(Tf-T,) 


(11) 

(12) 

(13) 
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The  minimal  Vbhp  depends  only  on  required  oxygen  molar  flow  rate  and  on  Tf  (or  Po).  If  N M i  defined 
by  (10)  is  very  low  or  less  than  0  the  BHP  volumetric  rate  should  be  higher  than  defined  from  water 
vapour  management  (13).  We  will  suppose  in  this  case  that  the  BHP  residence  time  should  be  short 
enough  to  support  NHf>Nu/2  and  good  chlorine  absorption  conditions.  Hence  the  maximum  BHP 
residence  time  is  in  this  case 

^BHP  =  (Nh/4)(H/M0)  (14) 

For  RSSOG  parameters  in  Tables  6  and  6a,  we  obtain  the  required  BHP  residence  time  in  RSSOG, 
consumption  of  BHP  alkaline  molarity,  the  axial  length,  axial  velocity  and  volumetric  rate  of  the  BHP. 
Hence  the  final  designing  parameters  and  expected  output  parameters  are  presented  in  Table  7. 


Tab.  7 

The  preliminary  designing  parameters  and  expected  output  parameters  for  M(O2)=30  mmole/s 

(explanation  in  Fig. 6) 


Input  parameters  (T=253K,  Po=100  torr,  Table  6) 

Output  parameters  (Tf=315K,  r|w<0.15) 

Item 

n 

s4 

H 

mm 

do 

(am 

Ug 

m/s 

Nh 

mole/L 

Mc 

mmole/cm2/s 

tBHP,S 

z 

cm 

w 

cm/s 

Vbhp 

L/s 

NHf 

mole/L 

1 

60 

3 

200 

130 

6 

~1 

0.5 

1.6 

3.2 

0.018 

2.7 

2 

100 

1.5 

100 

130 

6 

~1 

0.25 

1.6 

6.4 

0.018 

2.7 

In  Items  3-4  the  BHP  volumetric  flow  rate  should  be  higher  than  defined  from  (13)  or  tui n>  should  be 
defined  from  (14)  to  keep  NHf>NH/2  mole/litre 

3 

80 

6 

300 

130 

6 

~3 

0.3 

0.53 

1.8 

0.02 

3 

4 

80 

1 

50 

100 

3 

~0.5 

0.15 

3.2 

21.3 

0.04 

1.5 

5 

100 

1 

50 

50 

2 

-0.5 

0.1 

3.2 

32 

0.06 

1 

Input  parameters  for  the  moderate  pressure  RSSOG 
(T=253K,  P0=20  torr,  Table  6a) 

Output  parameters  (Tf=284K,  r|w  <0. 15). 

6 

80 

5 

500 

200 

6 

-3 

0.14 

0.53 

3.8 

0.036 

4.3 

7 

80 

4 

300 

100 

3 

-2 

0.1 

0.8 

8 

0.06 

2 

The  efficiency  of  the  BHP  consumption  at  high  Po  is  close  to  50%.  At  moderate  pressure  Po=20  torr 
the  consumption  of  the  BHP  is  lower.  The  most  suitable  input  parameters  are  in  Item  3  at  Po=100  torr 
and  Item  6  for  Pq=20  torr. 


Fig.  6  Designing  parameters  of  RSSOG 
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3.  Principle  constructions  of  the  small-scale  RSSOG  for  production  of 

30  mmole/s  of  oxygen  flow  (close  coupling  of  the  generator  to  nozzle  plenum  to 
minimize  02(1A)  loss  and  disengagement  of  gas/iiquid) 

The  volumetric  rate  of  the  gas  over  BHP  bubbler  layer  is 


v„,=”  (.5) 

'  0 

Rg=8.31  J/mole/K  is  the  universal  gas  constant.  The  gas  transport  time  from  the  BHP  bubbler  layer  to 
the  throat  of  the  nozzle  bank  is  equal  to 


_  Vo  _  VqPq 

Vgas  M(02)RgT  ’ 


(16) 


Vo  is  the  volume  of  the  transport  duct  between  the  BHP  surface  and  the  nozzle  throat. 

The  so-called  quenching  parameter  Px  due  to  reactions  02(1A)+02(iA)— +C>2(32)  is  equal  to 


_  PqVq  VqPq 
Vgas  M(02)RgT 


(17) 


g 


Let’s  require  (Pottmns)  <Xmax“0.1  torrxs.  Hence  the  maximum  volume  of  the  transport  duck 

(18) 

For  the  input  parameters  presented  in  Table  7,  the  maximum  volumes  of  the  transport  duct  between 
the  bubbler  BHP  layer  and  the  nozzle  throat  are  presented  in  Table  8. 


V  <  Xmax  M<02>r  t 

P  P 

r0  r0 


Tab.  8 

(Items  1-7  are  the  same  like  in  Tables  6  and  6a) 


Output 
pressure  P0, 
torr 

Axial  length 
of  RSSOG  Z, 
cm 

Chlorine  flux  through  bubbler 
surface  Mc 
mmole/cm2/s 

Required  maximum  volume  of  transport  duct 
from  BHP  surface  to  nozzle  throat  Vo,  cm3 

1 

100 

1.6 

1 

4.8 

2 

100 

1.6 

1 

4.8 

3 

100 

0.53 

3 

4.8 

4 

100 

3.2 

0.5 

4.8 

5 

100 

3.2 

0.5 

4.8 

6 

20 

0.53 

3 

120 

7 

20 

0.8 

2 

120 

Table  8  shows  that  the  required  transport  duct  volume  is  too  small  for  Items  1-5.  It  is  quite 
difficult  to  imagine  an  assembling  of  RSSOG  with  R=3  cm  and  Z«1  cm,  and  with  the  nozzle 
bank  keeping  the  plenum  pressure  Po=100  torr  and  keeping  Ch^A)  yield  (Items  1-5). 
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Much  more  promising  is  the  scheme  of  COIL  with  output  oxygen  pressure  Po=20  torr  (Items  6- 
7).  Also  the  COIL  schemes  with  RSSOG  parameters  according  to  Items  (1-5)  with  output  pressure 
Po=100  torr  can  be  realized  using  oxygen  pressure  discharge  from  the  output  pressure  Po=100  torr  to 
the  plenum  pressure  Ppien=20  torr.  The  discharge  of  the  pressure  can  be  realized  immediately  after 
oxygen  outlet  from  BHP.  In  this  case  the  required  maximum  transport  volume  is  equal  to  120  cm1 
(like  for  items  6-7)  and  large  enough  for  the  arrangement  of  the  nozzle  bank. 


Scheme  1 

The  principle  scheme  of  RSSOG  with  output  and  plenum  oxygen  pressure  20  torr  is  presented 
in  Fig.7.  The  RSSOG  can  be  scaled  in  Z  and  R  directions.  The  appropriate  changes  of  parameters  w, 
Vbhp  should  be  made  according  to  formulas  ( 1 1)-(  14).  For  example,  increasing  Z  twice  forces  an 
increase  in  the  axial  velocity  w  and  the  BHP  volumetric  rate  twice.  The  principle  assembling  of  three 
or  several  RSSOG  is  presented  in  Fig.8. 


A 


Oxygen 

outlet 

30  mmole/s 


Chlorine  inlet 


J  Nonrotating 
parts 


Rotating 

parts 


I  I  BHP 
line 


|  |  Chlorine 

line 


Fig.  7  Principle  scheme  of  RSSOG 


17 


Fig.  8  Principle  of  several  RSSOG  assembling 


Let’s  present  some  technical  characteristics  of  the  assembling  of  several  RSSOG.  The  oxygen 
plenum  flux  (in  front  of  the  nozzle  bank)  is  Mo2=nMc27iR/L,  where  n  is  the  number  of  RSSOG,  L  is 
the  length  of  assembling  on  n  RSSOG.  For  Mc=3  (mmole/s)/cnr  the  oxygen  plenum  flux  of  the  order 
of  7  (mmole/s)/cm2  can  be  generated  by  this  set  of  small  scale  RSSOG.  The  ratio  of  the  total  oxygen 
molar  flow  rate  to  the  BHP  volumetric  rate  of  the  order  of  1.5  (mole/s)/(litre/s)  for  Item  3  and  0.8 
(mole/s)/(litre/s)  for  Item  6  can  be  achieved. 
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Scheme  2 

The  principle  scheme  with  output  pressure  100  torr  and  plenum  pressure  20  torr  is  presented  in 
Fig.  9  and  is  similar  to  the  scheme  in  Fig.  8  with  the  difference  that  a  perforated  wall  is  installed  over 
the  BHP  surface  to  discharge  the  oxygen  pressure  from  100  torr  to  20  torr. 


Chlorine  inlet 

Nonrotating 
parts 


Oxygen 

outlet  30  mmole/s 
20  torr 


Fig.  9  Principle  scheme  of  RSSOG  with  oxygen  pressure  discharge 
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Scheme  3 

The  assembling  of  several  RSSOG  in  vertical  Z  direction  is  presented  in  Fig.  10.  In  this 
scheme,  the  length  Z  is  of  the  order  of  several  cm  and  axial  velocity  w  of  the  order  of  several  ten 
centimetres  per  second.  The  worked-out  BHP  drains  to  the  common  collector. 


BHP  inlet 


Fig.  10  Principle  scheme  of  RSSOG  with  vertical  scaling 
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4.  Conclusions 


The  conception  of  the  centrifugal  bubble  SOG  is  based  on  understanding  that  the  mass  transfer 
phenomena  for  these  conditions  are  much  faster  than  in  any  another  type  mass-transfer  apparatus. 
Hence  very  efficient  production  of  the  oxygen  flow  with  high  Oil* A)  yield  and  chlorine  utilization  is 
possible.  The  estimation  of  the  gas-liquid  hydrodynamic  of  centrifugal  bubble  SOG  is  based  on  the 
literature  data  reviewed  in  Report  001  of  this  grant.  The  estimation  of  chlorine  utilization  and  OjC'A) 
production  is  based  on  the  modern  understanding  of  the  mass-transfer  phenomena  in  the  chlorine-BHP 
system.  Nevertheless,  we  understand  that  all  estimations  are  very  rough  because  the  real 
hydrodynamics  and  the  mass-transfer  processes  in  the  centrifugal  bubbler  SOG  are  much  more 
complex.  Very  simplified  model  of  the  processes  described  above  gives  a  very  approximate 
parameters  of  RSSOG  to  achieve  required  output  parameters  and  suitable  for  COIL  application. 
According  to  these  estimations,  the  oxygen  flux  Mo2~3  (mmole/s)  per  1  cm2  of  the  bubbler  surface  can 
be  produced  at  the  output  oxygen  pressures  20+100  torr.  A  possible  production  of  higher  oxygen  flux 
(with  high  02('A)  yield,  high  utilization  and  pressure  >20  torr)  is  a  question.  For  Mo2<3  mmole/s  and 
high  Po,  the  transport  volume  duct  between  RSSOG  and  the  nozzle  ban  (NB)  should  be  very  small  to 
save  02(1A),  so  it  is  quite  difficult  to  invent  efficient  RSSOG  +  nozzle  bank  system.  The  discharge  of 
oxygen  plenum  pressure  to  20  torr  would  make  possible  to  solve  the  problem  of  oxygen  transport  from 
the  RSSOG  to  the  NB.  Using  of  diluted  BHP  with  low  KOH  (NaOH,  LiOH)  concentration 
substantially  decreases  the  efficiency  of  chlorine  utilization  and  very  hardly  designed  parameters  of 
RSSOR  (a  small  nozzle  diameter,  a  high  centrifugal  acceleration)  are  required.  Several  schemes  of 
RSSOG  and  scaling  procedures  are  proposed  in  this  Report.  The  RSSOG  can  work  in  the  single  pass 
burn  down  with  efficient  BHP  consumption.  The  predicted  ratio  of  chlorine  molar  flow  rate  to  BHP 
volumetric  rate  is  of  the  order  of  l(mole/s)/(litre/s)  in  the  single  pass  burn  down  mode.  This  ratio  is 
much  higher  than  for  another  SOG  types.  The  scaling  of  RSSOG  is  evident  and  clear.  A  plenum 
oxygen  flux  in  front  of  the  large  scale  NB  of  the  order  of  7  mmole/cm2/s  can  be  produced  by 
assembling  several  RSSOG. 
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B.Concept(s)  investigated  in  the  COIL  Laboratory  in  Prague 


Two  concepts  of  the  advanced  generator  of  singlet  oxygen  were  proposed  employing  a  centrifugal 
force  for  gas/liquid  disengagement:  the  Centrifugal  Film  SOG  (CentFilmSOG)  and  Centrifugal  Spray 
SOG  (CentSpraySOG).  They  were  described  in  the  interim  Report  001  submitted  to  the  EOARD  on 
31  July  2004.  Theoretical  estimations  of  both  concepts  showed  that  the  CentSpraySOG  would  better 
fulfill  requirements  given  in  Objects  and  Goals  of  this  investigation,  and  therefore,  our  recent  research 
has  been  focused  on  more  detailed  computational  modeling  of  parameters  of  this  generator,  and 
designing  devices  for  its  experimental  investigation. 

Centrifugal  Spray  Singlet  Oxygen  Generator 
(CentSpraySOG) 

This  concept  is  based  on  principle  of  the  efficient  generation  of  singlet  oxygen  in  mixture  of 
chlorine  and  BFIP  spray  (mist),  from  which  liquid  droplets  are  efficiently  separated  by  a  centrifugal 
force  within  very  short  time.  The  mist  of  small-diameter  droplets  of  tens  microns  provides  a  large 
specific  surface  area  and  a  short  diffusion  path  of  PIO2  ions  from  the  bulk  BFIP,  resulting  in 
substantially  better  BFIP  utilization  than  in  any  other  types  of  chemical  generators  (jet-type,  disk- 
type).  Calculations  performed  at  the  USAF  by  Blauer,  et  al.  [1]  on  the  aerosol  generator  showed  that 
the  generator  operating  at  moderate  CFl'A)  pressures  (-20  Torr)  and  with  a  high  CF^A)  yield  (-70%) 
may  be  designed  with  very  restrictive  dimensions  to  limit  a  time  of  the  two-phase  flow  below  a  few 
milliseconds.  Our  new  proposed  concept  of  CentSpraySOG  employing  a  combination  of  the  spray 
SOG  and  centrifugal  separator  of  both  phases  may  overcome  this  problem. 

1.  Principle  design  and  operation  of  CentSpraySOG 

The  proposed  concept  is  schematically  shown  in  Fig.  1.  Operation  of  the  device  can  be  simply 
described  by  the  following  processes: 

•  An  atomizing  nozzle  (Fig.  2)  produces  a  spray  (mist)  of  BFIP  droplets  by  means  of 
pressurized  Cf  /Hc  (N2)  gas  mixture. 

•  Droplets  of  partially  exhausted  BFIP  are  separated  from  g/1  mixture  by  a  centrifugal  force 
of  the  rotating  blade  separator  within  a  few  milliseconds  (Fig.  3). 
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A  liquid-free  gas  (containing  Oof' A),  residual  CI2,  buffer  gas  and  residual  water  vapor) 
leaves  the  generator  through  the  exit  coupled  to  the  laser  nozzle  throat,  while  the  liquid  is 
collected  in  the  bottom  tank. 


Fig.  1  Scheme  of  principle  of  CentSpraySOG  device 


Fig.  2  Scheme  of  two-phase  spray  nozzle 


Fig.  3  Scheme  of  centrifugal  separator 


2.  Estimation  of  main  characteristics  of  the  spray  SOG 

Parametric  estimation  of  the  spray  generator  was  performed  for  conditions  listed  in  Table  1. 
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Tab.  1 

Main  parameters  of  the  reaction  system 


Characteristics 

Labeling 

Magnitude 

Gas  entering  the  generator 

Cl2  +  He 

1  :  4 

Total  generator  pressure 

Ptot 

40  kPa  (300  torr) 

Gas  temperature 

T 

260  K 

BHP  droplet  radius 

a 

10  pm 

Diffusion  coefficient  of  chlorine  in  gas 

Dcg 

3.3  x  10'4  nr  s'1 

Diffusion  coefficient  of  chlorine  in  BHP 

Dcl 

5  x  10'10  nf  s'1 

Diffusion  coefficient  of  HOC  in  BHP 

Dhl 

5  x  10'1IJ  nf  s'1 

Diffusion  coefficient  of  oxygen  in  BHP 

Dal 

5  x  10'10  nf  s'1 

Initial  HOC  concentration  in  BHP 

ChLO 

5  kmol  m-3 

Initial  Cl2  concentration  in  gas 

CCGO 

3.7xl0'3  kmol  m"3 

Rate  constant  of  reaction  Cl2  +  HOC 

kp 

5x10s  nf  kmol  's"1 

Striking  Cl2  probability  at  BHP  surface 

Yc 

0.0075 

Striking  02('Ag)  probability  at  BHP  surface 

Ya 

0.003 

Henry  constant  of  chlorine 

Hc 

1 

02('Ag)  lifetime  in  BHP 

Ta 

2  x  10'6  s 

Rate  constant  of  02('Ag)  dimolar  reactions  in  (g) 

£qg 

3.7xl04m3kmoT1s"1 

The  calculations  of  characteristic  times  of  partial  steps  in  the  whole  process,  i.e.,  the  characteristic 
time  of  HO2"  depletion  in  droplet  controlled  by  interfacial  crossing  or  reactive  diffusion  in  liquid,  the 
characteristic  time  of  Ch  depletion,  the  estimation  of  CB^A)  detachment  yield  and  02(*A)  loss  in  the 
gas  phase,  were  presented  in  detail  in  the  interim  report  of  this  grant.  Below  are  summarized  therefore 
only  conclusions  from  these  estimations. 

•  The  process  of  Ch^Ag)  generation  in  the  proposed  conditions  is  controlled  first  by  chlorine 
absorption  (interfacial  crossing)  (t  <  0.3  ms). 

•  After  this  first  period,  the  process  is  controlled  by  HCB  diffusion  from  the  droplet  interior. 

•  Chlorine  is  nearly  exhausted  within  3  ms. 

•  ~  40  %  of  HO2  ions  of  BHP  is  consumed  within  this  time. 

•  The  CTf1  Ag)  yield  at  the  generator  exit  decreases  with  time  from  0.73  (after  1  ms)  to  0.50 
(after  4  ms). 
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3.  Detailed  theoretical  description  of  CentSpraySOG 


Tab.  2 

List  of  main  symbols 


a 

Radius  of  droplet 

b 

Radius  of  gas  volume  surrounding  the  droplet 

Cif  (; i=C,HA\f=G,L ) 

Molar  concentration  of  HCL’,  CI2,  and  Cff'A)  in  gas  and  liquid  phase 

CisG  >  CisL 

Molar  concentration  at  the  g/1  interface  on  the  gas  or  liquid  side 

CHL0  ’  CCG  0 

Initial  (bulk)  concentration  of  HO2"  and  CL 

CCG  ’  CAG 

Averaged  concentration  of  CL  and  CL^A)  in  the  gas  phase 

Djf  ( i=C,H,A;f=G.L ) 

Diffusion  coefficient  of  species  in  gas  or  liquid  phase 

t 

Time  since  droplets  escaped  from  the  nozzle 

Dimensionless  Henry  constant  of  species  i 

a,  =  l/4vtj/j 

Mass  transfer  coefficient  of  species  i  across  g/1  interface 

Yi 

Stricking  probability  of  species  i 

Vt.i 

Thermal  velocity  of  species  i 

r 

Radial  coordinate  from  the  center  of  droplet 

Ui 

Utilization 

A  volumetric  fraction  of  the  liquid  in  the  aerosol  can  be  expressed  as 

f  .  .  V1 


X  = 


1  +  -- 
2  c, 


CGOj 


(1) 


where  u  =  U H  tU c.  If  the  process  mns  for  an  infinitely  long  time  the  ultimate  concentrations  for  the 
case  of  initial  stoichiometry,  i.e. 

2  Vg  ccgo  =  Vl  Chlo  (2) 


for  u  <  1  would  be 


CHLoo  —  (1  U)CHL0  ’  CCGu>  —  0  > 

(3) 

and  for  u  >  1 :  cHLa3  =  0  , 

CCG«>  =  (1  —  1  /  U)CCG 0  ' 

(4) 

3.1.  Mass  transfer  in  liquid  phase 

A  temporal  profile  of  the  HCL"  concentration  in  the  liquid  can  be  considered  for  two  limiting 

cases. 

3.1.1.  Initial  (film)  stage 

In  the  initial  stage,  the  diffusion-reaction  kinetics  is  limited  by  molecular  transport  at  the 
interface.  It  can  be  assumed  that  the  flux  of  HCL  towards  the  phase  interface  is  constant  in  time.  The 
corresponding  solution  is  then  (adaptation  of  eq.  1  in  [2],  p.  242) 
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.  3/  5r2 -3a2  2 a2  sin (ra„/a)  DHa„  t\  ... 

Chl=chlo-Jh  ~  +  — - - - —L^r^ - exp - (5) 

a  10 DHa  uHr  an  sin  an  ^  a 

where  an ,  n  =  1,2, . . . ,  are  the  positive  roots  of 

tan  a  =  a,  or  a  cos  a  -sin  a  =  0 ,  (6) 

and  Jh  is  the  molar  flux  of  HO2". 

The  first  root  of  eq.  (6)  is  al  =  4.4934 .  Higher  roots  were  found  using  the  Newton-Raphson  method 
with  a'2>  =  an_x  +  n  .  The  iteration  of  more  general  form  of  this  equation  will  be  discussed  later. 

Using  dimensionless  variables  y  =  r/a  and  r  =  tDH  /  a2 ,  eq.  (5)  can  be  written  as 


chl=chlo\1-0h  +  5}  3  -  exp(-g„2r)  ], 

1  10  sin  an 


where 


The  surface  concentration  (y=l)  is 


D  c 
^h'-h  0 


'  1  “  1  ^ 

C/m  =cHL0  \-</>H  3t  +  --22_j—' exp(-a2r)  =  cHL0[\- ^H&Hs{r)\  (9) 

^  L  5  «„  J 

The  significance  of  the  sum  of  exponentials  decreases  sharply  with  growing  r.  For  r  =  0.1  it 
contributes  less  than  3%. 

The  mean  concentration  averaged  from  r  =  0  to  r  =  a  is 

ra  2 

cH,  4nr~dr  t.  ;  f 

77  _  Jo  /im 


If  4^  «  a  ,  the  influence  of  curvature  can  be  neglected  and  the  case  can  be  treated  as  diffusion  into 
a  semi-infinite  solid  (adapt,  eq.  7  in  [2],  p.  75): 


C  HT  C  rri  n 


h  . 


Dh  V  n 


ADHt)-2erf\ 2^t\  ’ 


where  x  =  a-r  and 


2  I*00  c2 

erfc  =  —j=  [  e~^  di;  . 
sin  Jx 


3.1.2.  Final  (depleted  surface)  stage 

In  the  final  stage,  the  diffusion-reaction  kinetics  is  limited  by  the  transport  of  HO2"  towards  the 
phase  interface.  We  assume  cHsL  «  cHL ,  which  allows  an  approximation  clhl  «  0 .  An  adaptation  of  eq. 


4  in  [2],  p.  233  gives 
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2fl^(-l)"  .  nnr 


HL  HLO 


7TT 


I 


-sin - expi 

n  a 


f  DHn2n2t 


The  average  concentration  at  time  t  is  given  as  (adapt,  of  eq.  8.  in  [2],  p.  234) 

6  V  1  ( 

~tL— exP 


r  DHn2n2t ^ 


"-HL  '-HLO  _ 

n  »=i  n 

The  flux  at  the  surface  (r  =  a)  is  given  as 

D  c 

Ut-HL  _  ^H'-'HLO 


Jh  =  ~Dh 


dx 


a 


2^  exp 


f  DHn2K2t 


(13) 


(14) 


(15) 


3.1.3.  Matching  the  initial  and  final  solution 

An  approximation  for  the  diffusion  kinetics  in  the  whole  domain  can  be  obtained  by  combining 
the  initial  and  final  solutions.  The  initial  and  final  solutions  are  matching  at  the  time  tj,  for  which  two 
conditions  are  satisfied:  (i)  the  fluxes  are  equal,  and  (ii)  the  mean  concentrations  are  equal.  In  order  to 
achieve  that  solution,  a  start  of  the  final  solution  is  postponed  by  a  time  interval  At,  (At  <  td) . 
Parameters  td  and  At  are  adjusted  to  satisfy  the  above  two  conditions. 

The  computation  is  done  as  follows.  We  assume  that  the  initial  flux  is  given,  and  then  the 
dimensionless  form  can  be  used 


,  _  J  H,inia 

TH  ,im  ~  ~p.  • 

*-'HC HLO 

The  mean  concentration  from  the  initial  solution  at  the  matching  point  is 


-  HL,ini,d  C  HLO 


and  the  dimensionless  flux  computed  from  the  final  solution  is  given  as 

_  JH,fin,da 


H ,  fin, d 


=  2^]exp 


D  r 

^H^HLO  n= 1 


-n2;:2  ~Y{td  -  At) 


(16) 


(17) 


(18) 


By  setting <j)H fmd  =<fiHJniid,  eq.  (17)  can  be  iterated  to  obtain (t d  -  At).  The  mean  concentration  at  the 
matching  point,  computed  from  the  initial  solution,  is  given  as 


6  ^  1 

exP 


-  n2 it2  — j  (td  -  At) 
a* 


K  n= 1  n 

The  right-hand  side  can  be  evaluated  because  (t  d  -  At)  is  already  known. 
By  setting  cH  Jm  d  =  cll  lm  d ,  we  then  have 


=  CH,0$fin(td  “A0-  (19) 
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(20) 


2  ^  QftSKt  Af) 
t,  =  a  - - - 

3  DhKm 

The  delay  time  is  then  obtained  as  At  =  td  -  (td  -At). 


3.2.  Mass  transport  at  g/l  interface 

3.2.1.  Transport  across  the  interface 

It  is  assumed  that  temperature  at  both  sides  of  the  phase  interface  is  the  same.  The  net  flux  of  the 
species  i  from  the  gas  phase  to  the  liquid  phase  is 

(21) 

3.2.2.  Diffusion-reaction  kinetics  of  CI2 

We  use  here  the  coordinate  x  =  r  -  a  .  The  diffusion-reaction  equation  is 


h  =  a,(cisG  -cIsG)  =  l/4yivIJ(cisG-HjcisL) . 


dcr 


—  D  d  CL  b  r  r  ~n 

—  L-'CL  2  h-p'-HL'~CL  u 


(22) 


dt  dxz 

Solution  of  this  equation  can  be  found  when  approximating  cHL(x)  « cHsL  (i.e.  the  concentration 
variation  of  HChis  much  smaller  then  of  CI2) 

x 


CCL  ~  CCsL  eXP 


XB 


where 


The  flux  of  Ch  is  given  by 


Xp  = 


Dr 


kpCHsL 


id  =  -dcl  =  -°CL  ^  exp  TT"  =  -Col  4DakpCHsL  exp 

dx  Xp  Xp  Xp 


and  at  the  surface 


icsL  ~  Dcl  —  cCsLyj DCLkPcHsL  . 


X 


Combining  the  equations  (26)  and  (21),  we  obtain 

aCCCsG 


Jcs  ~ ' 


1  +  acHcXp  /  Dcl 


(23) 

(24) 

(25) 

(26) 

(27) 


3.2.3.  Diffusion-reaction  kinetics  of  Ch^A) 

The  second  Fick’s  law  for  singlet  oxygen  has  the  form 


5c  A 


dt 


=  ZX 


D2a 


dx 


■  +  kpcHLcCL  kgpC^p , 


(28) 


In  the  steady  state,  the  left  side  of  this  equation  is  nearly  zero. 
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Using  (23),  setting cHL (x)  «  cHsL ,  and  the  condition  ofc^*  — >  oo)  =  0,  we  can  finally  obtain  by 
integration  (28) 

cAi  =ccslX^V{x!  Xp)  +  {c^L  -  cCsL%)  exp(x  /  XQL  ) ,  (29) 

where  XQL  =  ^ D M  / kQL  ,  (30) 


x  =  ^L-\xPixfJ-\\ . 
D  \  1 


The  corresponding  flux  of  (^('A)  is 


t\  ocA  ..  ,c csl%  -i"  r'  a  va  r'/-./.  -i" 

7m  =  -  -O',  exp  —  +  y  exp  — 

OX  A  p  A  p  ^ 


)• 


Especially  for  the  phase  interface  (x  =  0) ,  the  flux  is  then 

1  1 


JasL  ~  DClCC 


-D..C. 


XPl  +  XP/XQL  M  **  X. 


QL 


and  Jasl  =  -«A (ca,g  -  ■ 

By  combining  eqs.  (33)  and  (34),  we  obtain 

1 


(31) 

(32) 

(33) 

(34) 


D  c 

^ cr'Csi 


J  As 


D  c 

_ ^ALuAsG 

v  X„  1 +  XPIXQL  XqlHaj 

'  y  ^ 

1  +  ^ 

X, 


where 


QL  J 

y  _  ^AL 
(*AHA 


(35) 


(36) 


The  length  XA  can  be  interpreted  as  a  thickness  of  hypothetical  additional  layer  of  liquid,  through 
which  (^('A)  must  diffuse. 


Using  eq.  (35) 


JcsL 


J  As 


D  c 

_ ^AL^AsG 

l  +  XP/XQL  XqlHaj 

'  x  A 
1  +  7*4 

V  X&J 


(37) 


3.2.4.  Simplified  evolution  equations  for  CI2  and  Chf'A)  concentrations  in  gas 

Because  the  diffusion  of  CI2  and  Cbf'A)  in  the  gas  phase  is  fast  in  comparison  with  the  H02" 
diffusion  in  the  liquid,  it  is  possible  to  ignore  the  spatial  gradient  of  their  concentrations: 

C  ~  C  /-• 


~'CG  ^ CsG  ’ 

CAG  ~  CAsG  ■ 

Then  the  evolution  of  chlorine  concentration  is  given  as 


(38) 
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(39) 


\n(b''  - ai)d<:^  =  Am2  jCs 
3  at 

dcrr  1  31 

—jr=~r^Ja- 

at  a  l- A 


It  follows  from  the  stoichiometry  that 


V  A 

CCG  (0  =  CCG0  +  2  y  [ CHL  (0  —  CHL0  ]  =  CCG0  +  [C HL  (0  —  ]  •  (41) 

The  evolution  of  (^('A)  concentration  is  given  as 

-a3)^jf-  =  4xa2jM  =  ^-x(b} -a3)kQGc2AG,  (42) 

3  at  3 

^  =  '  ^,j„-kQGc\c.  (43) 

at  a  l- A 

For  constant  flux  jA  and  initial  condition  cAG(t  =  0)  =  0  the  equation  (43)  can  be  integrated: 

c\G(t)  =  ,  1  tanh  1  ,  (44) 

k-QG*  a  4 


where  the  characteristic  time  is  tA 
Using  eq.  (43)  and  (37),  we  obtain 


1(1  -A)a 
I  34/ A shoe 


The  last  equation  can  be  written  as 


dcAG  1  31  [l  +  Xp/XQL  XqlHa)  f  2 
dt  a  1- A  r  „  \  WaG- 


1  + 


'  —  dp  +  +  Qgq  > 


where 


JcsL 

1  31  1  +  Xpl  XgL 

a  1  -1  f  X  1 


DaLCAsG 


1  31  XgLHA 

a  1-1 f  x 
1  +  ^- 
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QqG  ~  k QGCAG  • 

(50) 

After  formal  integration. 

CQG  ~  Qp  QqL  +  QqG  ’ 

(51) 

where 

Q,  =  \'oq,(t')dt'  ■ 

(52) 

4.  Numerical  solution  of  generator  parameters 

Differential  equations  described  in  the  chapter  3  were  solved  numerically  for  general 
parameters  listed  in  Table  3. 


Tab.  3 


Characteristics 

Label 

Magnitude 

Diffusion  coefficient  of  chlorine  in  BHP 

Dcl 

5  x  10"10  m2  s'1 

Diffusion  coefficient  of  HO2-  in  BHP 

DHl 

5  x  10'1U  mz  s'1 

Diffusion  coefficient  of  oxygen  in  BHP 

Dal 

5  x  10'1U  m2  s"1 

Rate  constant  of  reaction  CI2  +  HCB” 

kp 

5x10s  nr'kmol'V1 

Sticking  probability  of  CB  at  BHP  surface 

Yc 

0.0075 

Sticking  probability  of  Cff'Ag)  at  BHP  surface 

Ya 

0.003 

Henry  constant  of  chlorine 

Hc 

1 

Henry  constant  of  oxygen 

Ho 

20 

CBl'Ag)  lifetime  in  BHP 

Ta 

2  x  10-6  s 

Rate  constant  of  Cff'Ag)  dimolar  reactions  in  (g) 

kA 

3.7xl04  mj  kmof's"1 

Two  examples  of  calculated  results  for  various  input  parameters  are  shown  below. 

Example  1 


Input  Data 


Characteristics 

Label 

Magnitude 

Gas  entering  the  generator 

Cl2  +  He 

1  :  4 

Total  generator  pressure 

Ptot 

40  kPa  (300  torr) 

Diffusion  coefficient  of  chlorine  in  gas 

Dec 

1.04  x  10'4  nr  s"1 

Gas  temperature 

T 

260  K 

BHP  droplet  radius 

a 

10  pm 

Initial  HCB”  concentration  in  BHP 

CHL0 

3  kmol  m-3 

Ratio  of  BHP/chlorine  utilization 

Uh/Uc 

0.5 
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Results 


Characteristics 

Label 

Magnitude 

Volume  fraction  of  liquid  phase 

X 

0.0049 

Radius  of  gas  sphere 

b 

59  pm 

Initial  flux  of  chlorine  per  unit  droplet  surface 

jcsl 

1.7  mol  m"2  s"1 

Initial  characteristic  depth  of  Ch^Ag)  production 

XP 

0.58  nm 

Characteristic  depth  of  Cff'Ag)  quenching  in  liquid 

XQL 

31.6  nm 

Initial  flux  of  ChC'Ag) 

jvsl 

1.6  mol  m"2  s"1 

Duration  of  initial  (surface-limited)  period 

td 

0.24  ms 

Delay  of  the  stage  limited  by  H02  diffusion 

At 

0.13  ms 

Time  of  maximum  02(,Ag)  partial  pressure 

tmax 

6.2  ms 

Maximum  partial  pressure  of  02('Ag) 

PAmax 

5.5  kPa  (41  torr) 

Chlorine  utilization 

Uc 

100% 

BHP  utilization 

UH 

50% 

Fig.  4  illustrates  a  fast  HOj  depletion  from  the  BHP  surface.  In  this  period,  the  surface  process  plays 
a  controlling  role.  Afterwards  the  overall  process  is  controlled  by  HO2-  diffusion.  The  end  of 
decreasing  in  HO2-  mean  concentration  (at  t  =  6  ms)  corresponds  to  exhausting  time  of  chlorine  in  gas 
phase.  In  Fig.  5,  the  first  part  of  the  blue  curve  documents  an  accumulation  of  Chl'A)  during  the 
production  period,  and  the  declining  part  reflects  mainly  CB^A)  quenching  in  the  gas  phase.  It  means 
that  singlet  oxygen  should  exit  the  generator  in  a  time  of  the  maximum  CTI'A)  concentration,  i.e.  about 
6.2  ms  after  droplet  formation. 


Fig.  4  Time  dependence  of  H02‘  depletion 


Cl2  concentration 
_  02(A)  concentration 


3 


Fig.  5  Time  evolution  of  Cf  and  CBl'A)  concentration 
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Example  2 


This  case  was  calculated  for  a  higher  total  pressure  (100  kPa  (750  torr)),  lower  BHP  concentration, 
and  smaller  droplets  (d  =  2a  =  10  pm).  The  results  are  in  Figs.  6  and  7. 

Input  data 


Characteristics 

Label 

Magnitude 

Gas  entering  the  generator 

Cl2  +  He 

1  :  4 

Total  generator  pressure 

Ptot 

100  kPa  (750  torr) 

Diffusion  coefficient  of  chlorine  in  gas 

Dec 

0.33  x  10'4  m2  s'1 

Gas  temperature 

T 

260  K 

BHP  droplet  radius 

a 

5  pm 

Initial  H02~  concentration  in  BHP 

Chlo 

2  kmol  m-3 

Ratio  of  BHP/chlorine  utilization 

Uh/Uc 

0.5 

Results 


Characteristics 

Label 

Magnitude 

Volume  fraction  of  liquid  phase 

X 

0.018 

Radius  of  gas  sphere 

b 

19  pm 

Initial  flux  of  chlorine  per  unit  surface  of  droplets 

jcsl 

4  mol  m2  s"1 

Initial  characteristic  depth  of  02(*Ag)  production 

xP 

0.7  nm 

Characteristic  depth  of  02(’Ag)  quenching  in  liquid 

XQL 

31.6  nm 

Initial  flux  of  02('Ag) 

JasI 

3.8  mol  m"~  s"1 

Duration  of  initial  (surface-limited)  period 

td 

0.02  ms 

Delay  of  the  stage  limited  by  H02~  diffusion 

At 

0.01  ms 

Time  of  maximum  02(*Ag)  partial  pressure 

tmax 

1.4  ms 

Oxygen  pressure  at  this  time 

P02 

20  kPa  (150  torr) 

Maximum  partial  pressure  of  02(lAg) 

PAmax 

15  kPa  (112  torr) 

Chlorine  utilization 

Uc 

97% 

BHP  utilization 

UH 

49% 
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Cl2  concentration 

0.009 [  -  Oj(A)  concentration 

0.008 


Fig.  6  Time  dependence  of  HO;"  depletion 


Fig.  7  Time  evolution  of  Cl;  and  CTf'A)  concentration 


A  comparing  these  results  with  results  of  the  Example  1  shows  that  in  this  case,  the  initial  flux  of 
Chf'Ag)  and  its  maximum  partial  pressure  are  significantly  higher,  and  all  characteristic  times  are 
substantially  shorter.  A  maximum  partial  pressure  of  0;( '  Ag)  of  15  kPa  (112  torr)  should  be  achieved 
1.4  ms  from  the  beginning  of  chlorine  contact  with  BHP  droplets.  At  this  time,  the  calculated  CTf'Ag) 
yield  is  75%. 


5.  Two-phase  nozzles  for  spray  formation 

General  considerations 

The  preliminary  calculations  showed  that  high  partial  pressures  of  singlet  oxygen  (>  5.5  kPa 
(40  Torr))  could  be  achieved  with  only  very  small  droplets  in  formed  BHP  spray  (diameter  <  20  pm). 
According  to  a  survey  of  commercial  nozzles  available  from  world  producers,  small-droplet  sprays  of 
fairly  viscous  liquid  (like  BHP  at  low  temperatures  (10  -  50  mPas)),  could  be  produced  by  nozzles 
atomizing  liquid  by  pressurized  gas.  The  air-blast  atomizers  have  the  advantage  of  thorough  g/1  mixing 
and  production  of  fine  sprays.  In  a  pre-filming  type  nozzle,  liquid  is  first  spread  out  in  a  thin 
continuous  sheet  and  then  subjected  to  atomizing  action  of  high-velocity  gas.  Two  separate  gas  flows 
allow  the  atomizing  gas  to  impact  on  both  sides  of  the  liquid  sheet.  The  nozzle  can  produce  either  a 
hollow  cone  spray  (for  a  shorter  contact  time)  or  a  full  cone  spray. 

The  crucial  factor  for  choosing  a  suitable  nozzle  for  the  proposed  CentSpraySOG  is  the 
ratio  of  gas/liquid  flow  rates  that  ensure  a  time  atomization.  To  obtain  the  highest  Cb  utilization, 
rather  low  g/1  ratio  is  needed,  but  to  ensure  a  fine  atomization,  the  gas  flow  rate  must  be  as  high  as 
possible.  To  exclude  a  precipitation  of  the  salt  from  BHP,  which  could  block  the  nozzle  orifices,  the 
nozzles  with  external  atomizing  (mixing  gas  and  liquid  phases)  are  preferred.  Further,  a  choice  of  the 
nozzle  depends  also  on  requirement  for  the  contact  time  of  both  phases,  which  is  determined  by  the 
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shape  of  the  spray  cone.  The  contact  time  is  shorter  in  a  hollow  cone  spray  than  in  a  full  cone  spray. 
The  manufacturers  give  in  their  prospectus  usually  only  data  on  the  liquid  and  gas  flow  rate  for 
different  input  pressures  of  liquid  and  gas  but  the  data  on  the  droplet  size  are  given  only  exceptionally 
(and  only  for  water  and  air  at  atmospheric  pressure  and  for  selected  flow  rates  only). 

From  the  performed  survey  of  available  nozzles,  we  consider  following  types  of  the  nozzles  for 
our  planned  application  in  the  investigation  of  CentSpraySOG: 

BETE  Fog  Nozzle,  Inc.  (50  Greenfield  Street,  Greenfield,  MA  01301) 

The  SAM-type  nozzles  with  external  mixing  were  designed  for  spraying  of  viscous  liquid  yielding 
very  fine  droplets  (no  data  on  droplets  diameter  are  given). 

Tab.  4 

Gas  (atomizing  +  fan)  flow  rate  to  liquid  flow  rate  ratio  for  SAM  nozzles 


Type 

V G(atom)  /  Vl, 

(V  G(atom)  +  VG(fan))  /  Vl, 

(flG(atom)  IlG(fan))  /  Vl? 

mmol  (g)/ml  (1) 

SAM  05-03 

0.55 

0.82 

13.6 

SAM  06-04 

0.27 

0.42 

6.98 

SAM  07-05 

0.166 

0.25 

4.15 

Dusen-Schlick  GmbH  (Hutstrasse  4,  96253  Untersiemau/Coburg,  Germany)  Nozzles  of  the 
series  938,  848,  0/28,  and  0/48  can  atomize  liquids  to  a  finer  degree  than  other  two-substance  nozzles; 
models  0/48  and  840  are  designed  for  double  air  atomization. 

Tab.  5 

Atomizing  gas  flow  rate  to  liquid  flow  rate  for  two-substance  nozzles  providing  finest  atomization 


Type 

VG, 

vl, 

nG/vL, 

Nm3/h 

1/min 

mmol  (g)/ml  (1) 

938,  0/28 

11.0 

0.9 

9.1 

938,  0/28 

19.0 

1.5 

9.4 

0/48 

28 

0.7 

30 

0/48 

74 

2.5 

22 

0/48 

88 

4 

16 

Water  droplet  diameter  of  3.5  pm  to  27  pm  is  given  for  nozzle  Model  0/48  at  water  flow  rate  of 
0.4  1/min,  and  airflow  rate  of  74  Nm3/h  with  a  spray  cone  angle  of  120°. 
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Spraying  Systems  CZ,  s.r.o.  (Kamenna  31,  639  20  Brno,  Czech  Rep.) 

Tab.  6 

Type  1/8J  JAU,  Internal  Mix.,  no  data  on  droplet  size 


Type 

VG, 

VL, 

nG/vL, 

Nm3/h 

1/min 

mmol  (g)/ml  (1) 

SUJ  22DF 

3.4 

0.52 

6.6 

SUJ  22DF 

6.2 

1.33 

3.45 

For  the  spray  SOG  operation  using  chlorine-helium  mixture  in  the  ratio  1 :4,  and  2  to  8  M  ITCF" 
in  BF1P,  the  ratio  nG/vL  of  2.5  to  10  mmol/ml  is  required  for  attaining  the  BFIP  utilization  of  0.5. 

6.  Model  calculations  for  chosen  nozzles 

6.1.  Results  of  parametric  estimations 

The  model  calculations  for  the  nozzles  SAM-05-03  to  07-05  with  external  g/1  mixing  and  two 
gas  inlets  (for  atomization  gas  and  gas  controlling  the  spray-cone  angle)  were  performed.  The 
examples  below  were  calculated  for  gas  and  liquid  flow  rates  that  assure  a  fine  liquid  atomization.  The 
pressure  was  calculated  from  the  molar  flow  rate  of  gas  and  pumping  speed  corresponding  to  either 
the  rotary  pump  capacity  (300  m7h),  or  this  capacity  lowered  by  2x,  3x,  4x,  and  5x,  respectively.  For 
these  conditions,  time  dependences  of  the  surface  and  mean  concentration  of  HCF’  ions  in  the  drop  and 
the  mean  concentration  of  chlorine  and  singlet  oxygen  in  the  gas  were  calculated.  They  are  shown  in 
Figs.  8  -  11.  The  relative  concentration  of  F102  ions  in  Figs.  8  -  10  is  related  to  the  initial  HCF 
concentration  in  the  bulk  of  BFIP,  the  relative  CF  and  CF^A)  concentrations  in  gas  are  related  to  the 
initial  CF  concentration  in  the  gas  bulk. 

Figs.  8  and  9  illustrate  time  dependences  calculated  for  the  nozzle  SAM  06-04  and  for  two 
different  droplet  diameters  in  the  spray.  It  can  be  seen  that  1TCF  ions  on  the  drop  surface  are  quickly 
exhausted  (within  1.5  ms).  The  decrease  in  the  mean  HO2"  concentration  in  liquid  finishes  after  the 
chlorine  exhausting.  Neglecting  the  effect  of  limited  chlorine  transport  rate  in  gas  causes  a  sharp 
change  in  the  course  of  CF  concentration  near  zero.  A  maximum  of  CFf'A)  concentration  was 
achieved  after  24  and  13  ms,  respectively.  The  CF('A)  yield  was  still  rather  high  at  this  time  (Y  =  0.64 
or  0.7,  respectively).  The  CF  exhausting  was  several  times  faster  in  the  finer  spray  (Fig.  9)  where  the 
specific  surface  area  of  liquid  was  4times  higher.  The  CF  and  BFIP  utilizations  were  also  higher  for 
the  spray  with  smaller  drops  (Uc  =  0.75  and  Uh  =  0.38). 
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Fig.  8  Time  dependence  of  relative  mean  (1)  and  surface  (2)  concentrations  of  H02 
ions  for  40  pm  drop,  mean  concentration  of  Cl2  (3)  and  O^’A)  (4)  in  gas, 
and  02('A)  yield  (5),  cHo2-  =  5  M,  vL  =  31.7  ml/s,  ncl2  =  40  mmol/s  Cl2, 
Cl2:He  =  1:4,  Qpump  =  0.016  m3/s,  Ptot  =  32  kPa  (240  torr) 


Fig.  9  Time  dependence  of  relative  mean  (1)  and  surface  (2)  concentrations  of  H02  ions  for 
20  pm  drops,  mean  concentration  of  Cl2  (3)  and  02('A)  (4)  in  gas,  and  02('A)  yield  (5); 
other  conditions  as  in  Fig.  8. 


Results  of  calculations  for  a  larger  nozzle  (SAM  07-05)  and  higher  chlorine  flow  rate  (and  so  a  higher 
pressure),  are  in  Fig.  10. 


37 


Fig.  10  Time  dependence  of  relative  mean  (1)  and  surface  (2)  concentrations  of  H02“  ions 

for  20  pm  drops,  the  relative  mean  concentration  of  Cl2  (3)  and  02(*A)  (4)  in  gas,  and  02('A) 
yield  (5),  nCi2  =  60  mmol/s  Cl2,  CkHe  =  1:4,  cHo2  =  5  kmol/m3,  vL  =  58.3  ml/s, 

Ptot  =  48  kPa  (360  torr) 


Fig.  11  Time  dependence  of  relative  mean  (1)  and  surface  (2)  concentrations  of  H02  ions 

for  20  pm  drops,  the  relative  mean  concentration  of  Cl2  (3)  and  02('A)  (4)  in  gas,  and 
O^'A)  yield  (5),  nCi2  =  60  mmol/s  Cl2,  Cl2:He  =  1:4,  cHo2  =  3  kmol/m1,  vL  =  58.3  ml/s, 
Ptot  =  48  kPa  (360  torr) 
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The  Ch  utilization  is  higher  in  this  case  (Uc  =  0.81)  than  in  Fig.  9,  which  is  caused  mainly  by  a  higher 
surface  area  of  liquid  due  to  the  lower  g/1  ratio  for  this  nozzle  (see  Table  4).  The  BHP  utilization  is 
slightly  lower  (Uh  =  0.34)  than  in  Fig.  9,  but  it  could  be  higher  for  less  concentrated  BFIP 
(3  M  HO2")  (Fig.  11).  The  chlorine  utilization  is  then  however  lower. 

It  should  be  mentioned  that  the  results  of  calculations  above  are  rather  tentative  because  of  the 
real  drop  size  of  the  BHP  produced  by  these  nozzles  is  approximate  only,  and  more  reliable 
information  about  this  parameter  must  be  obtained  experimentally.  Planned  measurements  on  the  drop 
size  will  be  performed  first  with  a  model  solution  of  the  same  viscosity  like  BHP  (NaOH  or  glycerin), 
and  then  with  BHP  at  the  temperature  used  for  the  COIL  operation. 

Pressure  dependence  of  the  spray  SOG  parameters 
Example  1 

This  calculation  was  performed  for  different  CL  flow  rates  and  one  BHP  flow  rate  using  the 
nozzle  SAM  06-04.  The  results  correspond  to  the  time  of  maximum  O 2 ( '  A )  partial  pressure  and  are 
shown  in  Figs.  12  a,  b,  c. 


Fig.  12  a 
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Figs.  12  a,  b,  c  Dependence  of  chlorine,  BHP  utilization,  and  02(1A)  yield  on  total  pressure; 

nozzle  SAM  06-04,  d  =  40  pm,  vL  =  3 1 .7  ml/s,  ncl2  =  30-60  mmol/s  Cl2, 
xCi2  =  0.2  (Cl2:He  =  1:4),  cH02-=  5  krnol/m3 


A  decrease  in  Uci2  with  increasing  pressure  can  be  explained  by  a  higher  CB^A)  loss  in  gas  at  higher 
pressure.  At  higher  pressure,  the  maximum  concentration  of  CB^A)  is  achieved  at  lower  degree  of  Cl2 
conversion  to  Ch^A).  Chlorine  utilization  Uc  declines  markedly  with  increasing  n(  i2  at  constant 
pumping  velocity.  This  decrease  is  much  stronger  than  it  would  be  correspond  to  the  pressure  effect. 
The  BHP  utilization  (Uh)  weakly  increases  with  increasing  Cl2  flow  rate.  The  ChC'A)  yield  decreases 
with  pressure  at  p<25  kPa  (-200  Torr).  A  negligible  effect  of  pressure  on  the  ChpA)  yield  observed  at 
pressures  above  25  kPa  (Fig.  12c)  is  difficult  to  explain  because  every  point  relates  to  a  maximum 
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02(‘A)  concentration,  not  to  either  constant  reaction  time  or  reactor  volume.  In  this  example,  the 
highest  partial  pressure  of  oxygen  was  4-5  kPa  (30  -  37  Torr)  at  medium  (^('A)  yield  of  0.68. 


Example  2 

This  calculation  was  performed  for  various  dilutions  of  chlorine  with  helium  with  constant 
chlorine  flow  rate  in  the  nozzle  SAM  06-04.  At  a  lower  molar  fraction  of  chlorine,  the  He  flow  rate  was 
higher  resulting  in  a  higher  total  pressure,  while  Uc  remained  constant.  This  is  due  to  neglecting  mass- 
transfer  constrains  in  the  gas  phase  applied  in  this  model.  The  results  are  presented  in  Figs.  13  a,  b. 


total  pressure,  Pa 


total  pressure,  Pa 


Fig.  13  a,  b  Dependence  of  chlorine  utilization,  and  CTf'A)  yield  on  total  pressure; 

nozzle  SAM  06-04,  d  =  40  pm,  vL  =  3 1 .7  ml/s,  nCi2  =  40  mmol/s  CI2, 
xCi2  =  0.1  -  0.25,  i.e.,  nHe=  120  -  360  mmol/s,  cHo2-=  5  kmol/m3 
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Example  3 

In  this  calculation,  the  effect  of  BHP  concentration  was  examined.  The  results  are  shown  in 

Figs.  14  a,  b. 


cBHP,  mol/I 


cBHP,  mol/l 

Fig.  14  a,  b.  Dependence  of  chlorine  utilization,  and  CpC'A)  yield  on  BHP  concentration 
at  two  different  pressures  ;  nozzle  SAM  06-04,  d  =  40  pm,  vL  =  31.7  ml/s, 
nCi2  =  40  mmol/s  CP, 

The  increasing  Uc  with  increasing  HCB  concentration  is  caused  by  the  effect  of  HO2  mass  transfer 
rate  controlling  the  rate  of  CTl' A)  formation  in  the  final  stage  (see  the  chap.  3).  The  CTf'A)  yield  is 
nearly  independent  of  the  HO2  concentration. 
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Example  4 

This  example  is  calculated  for  a  larger  nozzle  (SAM  07-05)  making  possible  higher  flow  rates 
of  BHP  at  the  same  gas  glow  rate,  and  so  a  higher  ratio  of  liquid/gas  volume.  The  results  are  in  Figs. 
15  a,  b. 
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total  pressure,  Pa 


Fig.  15  a,  b  Dependence  of  chlorine  utilization,  and  CEf'A)  yield  on  total  pressure; 

nozzle  SAM  07-05,  d  =  40  pm,  vL  =  58.3  ml/s,  nCi2  =  30-  60  mmol/s  Cb, 
xCi2  =  0.2,  cH02-=  5  kmol/m3 

In  comparison  with  Example  1,  a  higher  Uc  can  be  achieved  for  the  same  gas  flow  rate  and  pressure. 
This  can  be  explained  by  a  higher  number  of  droplets  and  their  larger  specific  surface.  The  O2  partial 
pressure  of  5  -  7.7  kPa  (38  -  58  Torr)  with  Ch^A)  yield  of  0.65  -  0.8  was  calculated. 
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Example  5 

This  calculation  shows  the  effect  of  droplet  size  in  the  spray  on  the  CI2  utilization  and  Ch^A)  yield 
(4times  higher  specific  surface  area  of  liquid  vs.  Example  1).  The  results  are  in  Figs.  16  a,  b,  c. 

Fig.  16c  illustrates  the  dependence  of  effective  reactor  volume  on  the  pressure  at  CL  molar  flow  rates 
that  could  be  sufficient  for  1  kW -class  COIL.  It  can  be  seen  that  the  reactor  volume  is  relatively  small 
(up  to  100  cm3)  at  pressures  of  24  -  48  kPa  (180  -  360  Torr).  The  cylindrical  reactor  of  5  cm  in 
diameter  could  be  then  only  5  cm  long.  The  calculated  O2  partial  pressure  would  be  4.8-  9  kPa  (36  - 
67  Torr)  at  the  02('A)  yield  exceeding  0.8. 


total  pressure,  Pa 


total  pressure,  Pa 
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Fig.  16  a,  b,  c  Dependence  of  chlorine  utilization,  02(*A)  yield,  and  reactor 
volume  on  total  pressure;  nozzle  SAM  07-05,  d  =  20  pm, 
vL  =  58.3  ml/s,  nCi2  =  30-  60  mmol/s  Cl2,  xCi2  =  0.2,  cH02-=  5  kmol/m3 


Limitation  of  the  present  model 

The  model  used  for  the  above  calculations  was  based  on  several  assumptions,  which  may  result 
in  differences  between  its  prediction  and  experiment.  They  are: 

1/  Infinite  mass  transport  rate  in  the  gas  phase: 

This  assumption  becomes  incorrect  when  the  chlorine  utilization  reaches  100%.  At  that  time,  the 
transport  limitation  in  gas  arises  and  the  real  utilization  would  be  slightly  lower  then  predicted. 

2/  Perfect  “co-flow”  of  gas  and  liquid  phases 

In  the  real  situation,  the  velocity  of  gas  and  liquid  droplets  can  be  different,  especially  during  the 
initial  stage.  Droplets  could  also  coagulate  to  some  extent.  Both  these  effects  can  result  in  a  surface 
renewal  and  so  enhancement  of  the  mass  transfer  towards  and  backwards  the  droplet.  This  effect 
causes  a  faster  CB  utilization  and  can  thus  extend  a  time  for  singlet  oxygen  quenching.  Quantitative 
evaluation  of  the  impact  of  surface  renewal  on  the  generator  performance  can  be  estimated  only 
experimentally. 

3/  The  rate  of  CBt'A)  quenching  in  the  liquid  phase 
The  above  calculations  included  the  value  of  CBC’A)  lifetime  in  BHP  of  2  x  10'6  s  given  in  the 
literature.  Blauer  at  al.  [1]  found  by  comparison  between  their  model  and  experiments  that  the  liquid 
phase  quenching  could  be  higher  by  two  orders  of  magnitude  than  corresponds  to  the  literature  value. 
Based  on  this  estimation,  the  CBC'A)  yields  evaluated  from  our  modeling  can  be  treated  as  the  upper 
limit  values. 
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4/  Uniform  droplet  size 

The  real  spray  will  have  a  certain  distribution  in  the  droplet  size.  The  results  on  CL  utilization  and 
yield  are  ‘'‘correct”  for  the  mean  droplet  radius  included  in  the  model  calculations,  but  the  droplet  size 
distribution  must  be  taken  into  account  when  the  model  counts  also  with  the  enthalpy  balance.  The 
present  model  does  not  predict  yet  the  heat  release  and  convection,  and  consequently  the  amount  of 
water  vapor  content  in  the  gas.  It  is  obvious  that  smaller  droplets  will  evaporate  faster  and  so  produce 
more  water  vapor. 

6.2.  Conclusions  from  the  modeling 

•  The  calculations  performed  for  chosen  nozzles  producing  the  BHP  spray  with  40-p.m  or  20- 
pm  droplets  show  that  Oof' A)  can  be  generated  by  the  CentSpraySOG  of  a  small  inner  volume  100 
cm3,  at  the  oxygen  pressure  of  4.8-  9  kPa  (36  -  67  Torr),  and  with  the  CL^A)  yield  exceeding  0.8.  The 
SOG  with  these  parameters  would  be  satisfactory  for  driving  a  1  kW-class  COIL.  The  above 
calculations  were  performed  for  rather  concentrated  BHP  solutions  (with  3-7  M  HCL)  because  the 
chosen  small  two-phase  nozzles  need  a  relatively  high  g/1  ratio. 

•  Even  higher  CL^A)  partial  pressure  (20  kPa  (150  torr))  could  be  achieved  employing  a  spray 
with  smaller  droplets,  but  production  of  such  sprays  could  be  difficult  due  to  the  high  BHP  viscosity 
at  low  temperature.  This  problem  could  be  solved  by  using  higher  flow  rates  of  more  diluted  BHP, 
which  would  be  advantageous  also  for  the  lower  temperature  increase  during  the  reaction  with 
chlorine.  A  needed  lowering  the  gas/liquid  ratio  can  be  obtained  easier  by  using  larger  nozzles. 
Temperature  of  the  two-phase  mixture  could  be  lowered  also  by  some  cooled  liquid  (e.g.,  a  cooled 
diluted  H2O2  or  some  non-freezing  and  non-reactive  solution),  which  would  be  sprayed  from  another 
nozzle  downstream  of  the  reaction  region.  After  this  cooling,  the  two-phase  mixture  would  enter  the 
centrifugal  separator. 

7.  Experimental  set-up  for  measuring  a  drop  size  and  drop  velocity 

7.1.  Methods  available  for  these  measurements 

Methods  employed  in  a  drop  size  measurement  may  be  grouped  into  three  categories: 
mechanical,  electrical  and  optical.  The  optical  methods  are  either  droplet  imaging  (photography, 
holography)  or  droplet  non-imaging  (e.g.  Phase  Doppler  Anemometry).  Some  of  them  are  briefly 
described  below. 
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High-speed  photography 

This  technique  is  rather  simple  utilizing  a  flashlight  or  laser  pulse  to  create  a  high  intensity  and  short- 
time  light  source.  Duration  of  the  flashlight  should  be  of  the  order  of  1  ps,  and  the  laser  pulse  of  the 
order  of  nanoseconds.  It  can  be  used  also  for  getting  the  information  about  drop  velocities.  This 
technique  is  not  very  convenient  because  of  time-consuming. 

Light-scattering  interferometry 

Drop  size  and  velocity  measurement  is  based  on  the  observation  of  light  scattered  by  drops  flying 
through  the  crossover  region  of  two  intersecting  laser  beams.  The  drop  velocity  is  calculated  as 

U  =  X  fD/(2  sin(0/2)) 

where  X  is  the  wavelength  of  laser  light,  t'n  the  Doppler  frequency,  and  0/2  laser  beam  intersection 
half-angle.  A  size  information  is  contained  in  the  relative  modulation  (visibility)  of  the  scattered 
signal.  Based  on  the  scalar  diffraction  theory,  the  relationship  between  drop  size  and  visibility  can  be 
expressed  as  Visibility  =  2  Ji(jiD/5)/7tD/5 

where  D  is  drop  diameter,  Ji  the  first-order  Bessel  function,  and  8  the  interference  fringe  spacing. 
Phase  Doppler  Particle  Analysis 

Theoretical  description  of  the  dual-beam  light  scattering  showed  that  the  spatial  frequency  of  the 
scattered  interference  fringe  pattern  at  off-axis  angle  is  inversely  related  to  a  drop  diameter.  The 
method  uses  pairs  of  detectors  at  fixed  spacing  in  the  image  of  interference  or  fringe  pattern.  The 
overall  measurement  size  range  is  from  1  pm  to  several  thousands  pm,  and  can  be  used  also  for  dense 
sprays. 

Malvern  Particle  Sizer 

The  device  uses  a  method  based  on  the  Freuenhofer  diffraction  of  parallel  beam  of  the  monochromatic 
light  by  a  moving  drop.  A  diffraction  pattern  is  formed,  in  which  part  of  the  light  is  diffracted 
depending  on  the  drop  size.  This  commercial  technique  is  very  effective,  relatively  simple  and 
reliable. 

The  instrumentation  of  the  optical  methods  is  rather  expensive  and  their  use  needs  some 
experience.  It  requires  calibration  with  highly  mono-dispersed  droplets  that  are  either  sprayed  from  a 
dilute  liquid  suspension  and  subsequently  dried  or  present  in  liquid  suspensions. 

We  plan  to  employ  initially  the  high-speed  photography  for  exploration  of  the  spray  produced 
by  different  nozzles.  In  the  next  investigation,  some  non-imagine  optical  method  would  be  probably 
necessary. 
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7.2.  Scheme  of  device  for  spray  parametric  study 

The  device  consists  of  a  cooled  liquid  reservoir  connected  by  a  valve  with  the  spray-nozzle. 
The  liquid  will  be  compressed  by  nitrogen  delivered  through  tube  in  the  upper  part.  The  nozzle  is  fed 
by  atomizing  gas  (N2  or  He)  producing  a  spray  in  the  low-pressure  chamber  evacuated  by  rotary  pump 
(300  m3/h).  Side  mirrors  of  the  chamber  for  the  droplet  size  investigation  by  the  diagnostic  beam  are 
purged  by  means  of  nitrogen  (or  He). 


Fig.  20.  Scheme  of  device  for  spray  studies 

1  -  pressurized  N2  (He) 

2  -  inlet  of  tested  liquid 

3  -  tank  for  tested  liquid 

4  -  nozzle 

5  -  atomizing  N2 

6  -  diagnostic  chamber 

7  -  diagnostic  window 

8  -  inlet  of  purging  N2 

9  -  diagnostic  beam 


8.  Design  of  gas/liquid  blade  separator  for  CentSpraySOG 

8.1.  Geometry  of  channels  between  next  blades  in  separator 

The  separator  is  formed  by  rotor  with  radial  spiral  blades.  The  centrifugal  acceleration  of  the 
blades  is  a  =  ©2  r  (53) 
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where  ra  is  angular  velocity  (rad/s)  and  r  is  radius.  The  angular  velocity  is  related  to  a  number  of 
revolutions  per  minute  (r.p.m.)  as 

a>  =  2n  r.p.m.  /  60  .  (54) 

The  centrifugal  acceleration  can  be  decomposed  into  a  component  tangential  to  the  blade 

at  =  a  cos  a  ,  (55) 

and  a  component  normal  to  the  blade 

a„  =  a  sin  a  .  (56) 

Here  a  is  angle  between  the  tangent  of  the  blade  and  the  radial  coordinate. 

It  is  proposed  that  the  blade  shape  is  spiral  with  polar  angle  <9,  being  a  function  of  the  radial 


coordinate  r 

9  =  9m  —  tan  a;n  -  arctan  (1/tan  ctin)  +  (p2  -  \)m  +  arctan  (  l/(/?2  -  1)1/2  )  (57) 

Symbols  9m,  oiin,  and  r;n  relate  to  respective  parameters  at  the  inner  radius  of  the  blade. 

Here  p  is  a  dimensionless  radial  coordinate 

P  =  r/q  ,  (58) 

where  q  =  rm  /  (1  +  tan2  a;n)1/2  (59) 

The  angle  between  the  blade  tangent  and  radial  coordinate  is 

a  =  arctan  ( p  -  1)1/2  .  (60) 

The  length  of  the  blade  contour  between  radii  rm  and  r  is 

L(r)  =  (r2  -  rin2)  /  2cj  .  (61) 

This  spiral  blade  has  two  remarkable  properties: 

i)  The  tangential  acceleration  is  constant  along  the  blade  surface 

at  =  ©2  q  ,  (62) 

which  ensures  a  driving  of  the  BHP  film  with  constant  force. 

ii)  The  height  h  of  the  channel  formed  by  adjacent  blades  is  constant 

h=  A9r  cos  a  =  A9q  (63) 

where  A<9is  a  difference  of  polar  angles  of  adjacent  blades  of  the  same  radius. 


8.2.  Hydrodynamics  of  flowing  film  along  blades 

The  initial  approximation  for  the  flow  of  liquid  along  the  blade  is  made  for  the  assumption  that 
most  liquid  is  separated  from  the  g/1  mixture  at  the  entrance  into  the  blade  space.  The  relation  for  the 
tangential  velocity,  v(x)  was  taken  from  [3]  by  replacing  the  gravitational  acceleration  g  with  the 
centrifugal  acceleration  a 

v(x)  =  [(p  a$  cos  a)/2|i]  [1  -  (x/c5)2]  =  [(p  co2  ?<f)/2|i]  [1  -  (x/c5)2] ,  (64) 
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where  x  is  the  coordinate  normal  to  the  blade  (ranging  from  0  at  the  liquid/gas  interface  to  8  at  the 
liquid/blade  interface),  p  is  dynamic  viscosity  of  the  liquid. 

Maximum  velocity  at  x  =  0  is 

Vmax  =  (p  »2  qS2)  /  2|i  .  (65) 

Average  velocity  is  vav  =  2/3  vmax  =  (p  ff>2  q<T)  /  3(0.  .  (66) 

Volumetric  flow  rate  per  one  blade  of  width  w  is 

V\  =  w  5  vav  =  (wp  co2  qc f)  /  3p  .  (67) 

Average  retention  time  rav  is 

rav  =  L(r)  /  vav  •  (68) 

Depending  on  the  Reynolds  number 

Re  =  48  vav  p  /  p  (69) 

three  regimes  of  the  flow  can  be  recognized: 

•  Re  <  20:  laminar  flow  with  negligible  ripping.  Here  the  analysis  is  fully  valid. 

•  20  <  Re  <  1500:  laminar  flow  with  pronounced  rippling.  Here  the  analysis  can  be  good  for 
approximation. 

•  Re  >  1500:  turbulent  flow.  The  analysis  is  not  valid. 

8.3.  Settling  (disengagement)  of  droplets 

Mean  displacement  due  to  the  Brownian  motion  is 

Xrms  =  (2  Dp  t)l/2  (70) 

where  the  diffusion  coefficient  of  the  particle  (drop)  is  given  as 

Dp=kTCc  /  3nr/dp  .  (71) 

Ccis  the  Cunningham  correction  factor,  expressing  the  velocity  slip  between  the  particle  and  gas.  The 
terminal  settling  velocity  of  a  particle  is 

VTS  =  PpdpQn  Cc  / 1 8  q .  (72) 

Average  distance  travelled  due  to  a  settling  is 

Xs=vrst  =  pPdp2anCct/l8  r/ .  (73) 

Efficiently  removed  are  particles  for  which 

Xs  ^  Xrms  ,  (24) 

having  a  diameter  dp>  [(kTrj)/(ppanCct)]1/5  (75) 
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8.4.  Geometry  of  liquid  separator  and  parameters  of  liquid  film 

Dimensions  of  the  rotor  of  the  centrifugal  blade  separator  are  given  in  Table  7 


Tab.  7 


Rotor  inner  diameter, 

din  =  50  mm 

Channels  number 

n  =  6 

Channel  width 

w  =  50  mm 

Channel  length 

1  =  100  mm 

Blade  revolutions  per  minute 

n  =  1000  -  16000  r.p.m. 

Parameters  of  liquid  film 

Some  examples  of  the  calculated  effect  of  BHP  flow  rate  on  the  film  thickness,  the  mean 
velocity  of  liquid,  and  the  BHP  retention  time  were  presented  in  Figs.  5  -  7  of  the  interim  Report  001. 
Table  8  gives  the  film  thickness  calculated  for  rotor  parameters  from  the  Table  7,  the  total  BHP  flow 
rate  of  58.3  ml/s  (assumed  for  the  nozzle  SAM  07-05  in  the  Examples  4,5  in  chap.  6.1.),  and  BHP 
viscosity  of  20  mPas. 

Tab.  8 

Thickness  of  BHP  film  on  the  separator  blades 


rpm,  min'1 

1000 

2000 

4000 

8  000 

16  000 

§,  mm 

0.4 

0.2 

0.12 

0.09 

0.06 

The  calculated  BHP  layer  in  Table  8  is  rather  thin  obviously  due  to  its  high  tangential  acceleration. 

Separation  of  droplets  from  the  spray 

Table  9  gives  the  calculated  minimum  diameter  of  droplets  that  should  be  separated  in  the 
designed  separator  in  dependence  on  rotation  speed.  A  height  of  the  free  space  in  the  channel  for  the 
gas  flow  is  1  mm  obtained  by  subtracting  the  film  thickness  from  the  total  channel  height.  The 
calculation  was  performed  for  the  total  gas  pressure  of  50  kPa  (375  Torr). 

Tab.  9 

Minimum  diameter  of  droplets  (dp,mi„,  pm)  separated  by  CentSpraySOG 


rpm,  min'1 

1000 

2000 

4000 

8  000 

16  000 

dp, min? 

2.8 

1.8 

1.1 

0.67 

0.33 

It  is  obvious  that  even  small  droplets  could  be  separated  very  efficiently  in  the  proposed  centrifugal 
separator  at  revolutions  above  4000  r.p.m. 
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8.5.  Calculation  of  gas  pressure  loss  in  the  channel 

•  The  mean  height  of  the  free  space  in  the  channel  for  gas  flow,  he  =  0.8  mm 

•  Total  volumetric  gas  flow  rate  (60  mmol  Cl2  (02)/s  and  240  mmol  Heprim/s)  at  temperature  260  K, 
and  total  pressure  50  kPa  (375  torr): 

VG,tot=  nRT/P  =  0.3  R  260/50000  =  1.3  x  10'2  m3/h 

•  Total  free  cross-section  for  gas  flow  in  6  channels  (w  =  50  mm,  hG  =  0.8  mm): 

AG,,ot  =  6  x  0.05  x  8xl0-4  m2  =  2.4  x  10~4  m2 

•  Gas  velocity  in  the  channel 

VG,1  =  Vc.tot  /AG,tot  =  54.2  m/s 

•  Reynolds  number  ReG,i  =  v  d  pG  /  pG  , 

where  gas  density  pG  =  P  Mav  /(RT)  =  50000  x  0.0096/(8.314  x  260)  =  0.222  kg/m3 

molecular  weight  Mav  =  (60  x  0.032  +  240  x  4xl0"3)/200  =  9.6  x  10"3  kg/mol 

gas  (O2+  4He)  viscosity  pG=  1.87xl0"5  Pa  s 

ReG,i  =  vdpG/|iG  =  54.2  x  08xl0-3  x  0.222  /  1.87xl0-5  =  515 
=>  laminar  flow 

•  Pressure  loss  in  a  channel  between  blades  (laminar  flow  in  a  rectangular  channel): 

APg,i  =  4  v  1  p/(w  hG  k) 

k  =  0.3  hG/w  =  6  x  10"3  (tabled  value  of  the  coefficient  0.3  is  valid  for  hG/w  =  0.02) 

The  estimated  pressure  1ossAPg,i  =  4  x  54.2  x  0.1  x  1.87xl0"5/  (0.05  x  0.8xl0'3  x  6xl0'3) 

=  1688  Pa  =  1.7  kPa 

is  acceptable,  being  about  3%  of  the  total  pressure  in  the  generator. 


Calculation  of  Ch^A)  loss  in  the  separator 

The  deactivation  of  Oil'A)  in  the  gas  phase  goes  through  energy  pooling 


02('a)  +  02(‘a)  -4  02(‘a)  +  02(3Z) 

kx  =  2.7xl0"17  cm3/s 

(76) 

and  self  quenching 

02('A)  +  02(*A)  -4  02(3Z)  +  02(3S) 

ky  =  1.7xl0"17  cm3/s 

(77) 

The  effect  of  02('A)  quenching  can  be  expressed  as 

Y"1  =  (Y0)4  +  ktot  x 

(78) 

where  Y  is  the  O^'A)  yield  at  the  beginning  and  Y0  the  yield  at  the  exit, 

ktot  kx  2  ky  , 

X  =  (P.ot)2n02  V/(ntotRT)2 


and 
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(79) 

(80) 


Data  for  the  separator  with  6  gas  channels:  Ptot=50  kPa,  no2=0.060  mol/s,  ntot=0.3  mol/s 
V=6  x  0.1  x  0.05  x  8xl0"4  =  2.4  x  10"5  m3  (volume), 

VG,tot=l-3xl0"2  m3/h  (volumetric  flow  rate),  and  T=260  K 
Gas  retention  time  in  a  channel  of  the  length  1=0.1  m 
Tg.i  =  1/vg,i  =  0.1  /  54.2  =  1 .84  ms 
The  factor  y  calculated  from  (80)  equals  to 

X  =  8.56  x  10'3  mol  s/m3  =  5.15  x  1015  molec  s/cm3, 
the  product  ktot  X  is  then 

ktotX  =  0.314 

and  the  CF^A)  yield  at  the  separator  exit  is 
Y  =  0.64. 

This  yield  value  is  acceptable,  although  the  Cb^A)  loss  should  be  still  reduced  (the  product  Pt  is  0.138 
torr  s).  Either  smaller  height  of  the  gas  channel  or  a  shorter  channel  length  could  reduce  the  CE^A) 
loss. 

9.  Design  concept  of  CentSpraySOG  device 


A  cross-section  of  the  rotor  of  separator  was  schematically  shown  in  Fig.  3.  The  perpendicular 
cross-section  of  this  separator  is  given  on  Fig.  21.  The  rotor  consists  of  two  profiled  wheels  connected 
(by  welding)  by  six  blades.  Channels  for  gas  and  BHP  flow  are  formed  between  these  blades  of  the 
curved  shape,  see  calculations  in  chap.  8.2.).  The  reactor  inner  space  is  separated  from  the  outer 
atmosphere  by  aerodynamic  seals  formed  by  a  narrow  labyrinth  gap  between  rotor  and  stator.  The 
rotor  is  connected  with  the  electromotor  by  a  shaft.  This  motor  on  one  side  and  a  bearing  on  the  other 
side  are  protected  by  purging  with  nitrogen  or  helium.  The  purging  gas,  which  flow  rate  is  much 
smaller  than  a  gas  flow  through  the  nozzle,  streams  through  the  labyrinth  seal  and  exits  into  the 
reaction  space.  The  shape  of  trajectory  of  droplets  is  affected  by  the  pressure  loss  in  the  channels 
between  rotor  blades.  The  higher  is  this  loss,  the  more  uniform  is  the  g/1  flow  along  the  blade  width. 
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Fig.  21  Scheme  of  CentSpraySOG 


Dimensions  of  the  CentSpraySOG  for  about  1  kW  COIL  are  given  in  the  Table  7.  The  inner 
diameter  and  length  of  the  reactor  rotor  can  be  significantly  increased  for  attaining  higher  flow  rates  of 
singlet  oxygen.  The  increase  in  the  rotor  length  is  limited  by  a  need  of  the  comparable  droplet  delay  in 
the  inner  space  for  various  trajectories. 
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Fig.  22  shows  schematically  a  possible  arrangement  of  several  modules  of  CentSpraySOG  with 
the  outlet  into  one  supersonic  laser  COIL  nozzle. 


Fig.  22  Multi-module  CentSpraySOG;  1-  spray  nozzle,  2-  separator  body 


10.  Conclusions  of  present  study  on  SprayCentSOG 

A  summary  of  the  results  follows  the  requirements  on  the  advanced  SOG  characteristics  (Objects  and 

goals,  p.  4). 

1.  “A  high  pressure  operation  with  ~  100  torr  O2  partial  pressure By  our  calculations,  this 
parameter  could  be  fulfilled  in  the  proposed  CentSpraySOG  operating  with  a  droplet  spray 
containing  drops  less  than  20  pm.  The  calculated  chlorine  utilization  should  be  above  0.95  and 
02('A)  yield  above  0.75  (see  Example  2  in  chap.  4,  Example  5  in  chap.  6.1). 

2.  “A  single  BHP  pass  burn  down The  BHP  utilization  of  0.34  -  0.38  can  be  attained  for  the 
5  M  solution  of  HO2  ions,  which  is  one  order  of  magnitude  higher  utilization  than  for  the  jet  SOG. 
Even  higher  utilization  could  be  attained  using  more  diluted  solutions,  which  needs  however  to 
employ  nozzles  for  higher  flow  rate  of  liquid  (i.e.,  a  lower  ratio  of  gas/liquid  flow  rate). 
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3. 


“ A  close  coupling  of  the  generator  to  the  nozzle  plenum  to  minimize  O2C A)  losses”: 

The  CentSpraySOG  is  designed  for  a  direct  coupling  to  the  laser  nozzle.  Our  calculations  were 
focused  also  on  estimation  of  the  optimal  reactor  and  separator  volume.  Fig.  16c  illustrates  that 
these  requirements  are  very  restrictive.  In  the  reactor  volume  of  ~100  cm3,  the  02('A)  yield 
exceeding  80%  could  be  achieved  at  chlorine  flow  rate  of  50  -  60  mmol/s  and  oxygen  pressure  of 
36  -  67  Torr.  The  calculated  overall  CB^A)  yield  at  the  exit  from  the  separator  including  the 
reactor  (100  cm3)  and  separator  (45  cm3)  was  64%. 

4.  “An  efficient  disengagement  of  gas/liquid  mixture By  our  calculations,  this  process  could  be 
efficiently  attained  by  a  centrifugal  blade  separator  designed  for  the  CentSpraySOG.  At 
revolutions  over  4000  r.p.m.,  droplets  less  than  ~1  pm  should  be  removed  from  the  gas  stream. 

5.  “ A  scalability  to  be  suitable  for  the  airborne  and  mobile  application One  module  of  the 
CentSpraySOG  can  be  scaled-up  to  increase  production  rate  of  CB^A)  by  simultaneous  increase  in 
the  gas  and  liquid  flow  rate,  the  reactor  volume  and  the  rate  of  generation  exhausting.  A  choice  of 
suitable  spray  nozzles  is  very  broad,  and  larger  two-phase  nozzles  can  even  deliver  a  higher 
liquid/gas  volume  ratio.  For  the  need  of  airborne  application,  more  modules  of  the  CentSpraySOG 
could  be  coupled  to  one  laser  nozzle  (see  Fig.  22). 
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